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AERONAUTIC SYMBOLS 
I. FUNDAMENTAL AND DERIVED UNITS 



Symbol 

Metric 

English 

Unit 

Abbrevia- 

tion 

Unit 

Abbrevia- 

tion 

Length.. __ 

Time 

Force __ 

l 

t 

F 

meter 

second. 

weight of 1 kilogram 

m 

s 

kg 

foot (or mile) 

second (or hour) 

weight of 1 pound 

ft (or mi) 
sec (or hr) 
lb 

Power. _ 

P 

V 

horsepower (metric) 


horsepower 

hp . 

mph 

fps 

Speed 

/kilometers per hour 

(meters per second 

kph 

mps 

miles per hour _ _ 

feet per second. 

2. GENERAL SYMBOLS 


Weight = 7 x 1 # 

Standard acceleration of gravity= 9.80665 ml s 2 
or 32.1740 ft/sec 2 


W 

Mass=— 

g 

Moment of inertia=m/c 2 . (Indicate axis of 
radius of gyration k by proper subscript.) 
Coefficient of viscosity 


v Kinematic viscosity 

p Density (mass per-unit volume) 

Standard density of dry air, 0.12497 kg-m~ 4 -s 2 at 15° C 
and 760 mm; or 0.002378 lb-ft“ 4 sec 2 
Specific weight of “standard” air, 1.2255 kg/m 3 or 
0.07651 lb/cu ft 


3. AERODYNAMIC SYMBOLS 


Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio, ^ 

True air speed 
Dynamic pressure, ~pV 2 


Lift, absolute coefficient C L =-i y 

go 


Drag, absolute coefficient C B = 


D 

qS 


Profile drag, absolute coefficient C'e 0 = ^’ 

Induced drag, absolute coefficient C D 

- i go 

p_, 

qS 


Parasite drag, absolute coefficient C Dv — 


C 

Cross-wind force, absolute coefficient 

gS 


I'M 

it 

Q 

p 

R 


a 

€ 

«0 

<xt 


7 


Angle of setting of wings (relative to thrust line) 
Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 


Reynolds number, 



where Z is a linear dimen- 


sion (e.g., for an airfoil of 1 .0 ft chord, 100 mph, 
standard pressure at 15° C, the corresponding 
Reynolds number is 935,400; or for an airfoil 
ol 1.0 m chord, 100 mps, the corresponding 
Reynolds number is 6,865,000) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) 

Flight-path angle 
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A METHOD FOR DETERMINING THE RATE OF HEAT TRANSFER 
FROM A WING OR STREAMLINE BODY 

By Charles W. Frick, Jr., and George B. McCullough 


SUMMARY 

A method jor calculating the rate of heat transfer from the 
surface of an airfoil or streamline body is presented. A com- 
parison with the results of an experimental investigation indi- 
cates that the accuracy of the method is good. 

This method may be used to calculate the heat supply necessary 
for heat de-icing or in ascertaining the heat loss from the fuselage 
of an aircraft operating at great altitude. 

To illustrate the method , the total rate of heat transfer from 
an airfoil is calculated and compared with the experimental 
result. 

INTRODUCTION 

The calculation of the rate of heat transfer from the surface 
of bodies of aerodynamic shape is a problem for which no 
explicit solution is known to exist, although considerable 
attention has been given to the heat transfer from a heated 
flat plate into both laminar and turbulent flow regimes. It 
has been shown in reference 1 that the transfer of heat from 
the surface of a hot plate into an air stream flowing over the 
the plate is primarily a boundary-layer problem. The exten- 
sion of heat-transfer theory to permit the calculation of heat 
flow for an airfoil or streamline body is therefore a matter of 
calculating the boundary-layer characteristics which may be 
done by the methods of reference 2. The present report 
makes use of the results of references 1 and 2 to extend lieat- 
transfer theory to wings and bodies for which pressure distri- 
butions may be calculated with accuracy. The derivation of 
the method is given in the appendix. 

A limited experimental investigation of the method was 
made to determine the accuracy of the method in calculating 
the local rate of heat transfer into both the laminar and the 
turbulent boundary layer of an airfoil, and also to obtain a 
check on the computed total rate of heat transfer for the 
wing. 

It is hoped that this method will facilitate a more accurate 
determination of the heat losses from wings in designing heat 
de-icing systems as well as from fuselages in the design of 
cabin-heating systems for aircraft operating at great altitudes. 

SYMBOLS 

The symbols used throughout this report and in the appen- 
dix are defined as follows: 
c wing chord 

c P specific heat at constant pressure 


h heat-transfer coefficient, Btu/sq ft, ° F, sec 

k heat conductivity 

L length of streamline body 

M l local Mach number, ratio of velocity just outside the 

boundary layer to local velocity of sound 
M 0 free-stream Mach number, ratio of velocity of free 
stream to the velocity of sound in the free stream 
q x local rate of heat transfer, Btu/sq ft, sec 
R c Reynolds number based on wing chord ( Vcjv ) 

R l Reynolds number based on body length ( VL/v ) 

r radius to surface of streamline body at any point 

along the axis 

r Y radius to surface of streamline body at point for 

which boundary layer is being computed 
s distance along the surface from the stagnation point 

T local temperature inside boundary layer, ° F, absolute 
T l local temperature outside boundary layer, ° F, 
absolute 

T 0 free-stream air temperature, ° F, absolute 
t local temperature inside boundary layer, 0 F 

t L local temperature outside the boundary layer, ° F 
t 0 free-stream air temperature, ° F 

t p surface temperature corrected for compressibility, ° F 

t e surface temperature, ° F 

t p — 1 0 heat-transfer temperature difference, ° F 
U local velocity just outside the boundary layer 
If local velocity just outside the boundary layer at 
point for which boundary layer is being computed 
u local velocity inside the boundary layer 

V free-stream velocity 

x distance along the chord from the leading edge for 

an airfoil, or along the axis for a streamline body 
y distance normal to surface 

a angle of attack, degrees 

P eddy heat conductivity 

7 ratio of specific heats 

8 l heat-transfer characteristic length for a laminar 
boundary layer 

8 t lieat-transfer characteristic length for a turbulent 
boundary layer 
e eddy viscosity 

f turbulent boundary-layer parameter (£=^pU 2 /t 0 ) 

0 momentum thickness of boundary layer 

[-Jo 


i 
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/x absolute viscosity 

v kinematic viscosity 

p air density 

<7 Prandtl number (c P iJ,/k) 

t 0 surface shear 

APPARATUS 

Tests to determine the rate of heat transfer from a wing 
were made in the Ames 7- by 10-foot wind tunnel employing 
a heated, two-dimensional airfoil model. 

DESCRIPTION OF MODEL 

The model was a 7-foot chord NACA 65,2-016 airfoil which 
completely spanned the 7-foot dimension of the wind tunnel 
(fig. 1). Ordinates are given in table I. 



Figure 1 — Heated airfoil model installed in the wind tunnel. 


* The forward portion containing heating lamps was of sheet- 
aluminum construction, and the trailing section of laminated 
wood. Construction details are shown in figure 2. Built 
into the model were 26 pressure orifices and 41 iron-constan- 
tan thermocouples. The thermocouples were in three 


chordwise bands; one at midspan, and two, 4 inches on either 
side extending to 41 percent chord. The external surface of 
the heated portion was bare metal and that of the unheated 
wood part was painted. The entire model was made aero- 
dynamically smooth by chordwise rubbings with successively 
finer grit sandpaper. All power, temperature, and pressure 
leads were brought out through the top of the model which 
was sealed closed to prevent the escape of convected heat. 

TABLE I 

ORDINATES OF 7-FOOT-CHORD NACA 65, 2-016 AIRFOIL 


X 

Percent Chord 

X 

Inches 

y 

Percent Chord 

V 

Inches 

0 

0 

0 

0 

.5 

.420 

1.202 

1.010 

.75 

.630 

1.423 

1. 195 

1.25 

1.050 

1.796 

1.509 

2.5 

2. 100 

2.057 

2. 106 

5.0 

4. 200 

3. 543 

2. 976 

7.5 

6. 300 

4.316 

3. 625 

10 

8. 400 

4. 954 

4.161 

15 

12.600 

5. 958 

5.005 

20 

16. 800 

6. 701 

5. 629 

25 

21. 000 

7. 252 

6. 092 

30 

25. 200 

7. 645 

6. 422 

35 

29. 400 

7. 892 

6. 629 

40 

33. 600 

7. 995 

6.716 

45 

37. 800 

7. 938 

6. 668 

50 

42. 000 

7. 672 

6. 444 

55 

46. 200 

7.184 

6. 035 

60 

50. 400 

6. 495 

5. 456 

65 

54. 600 

5. 647 

4. 743 

70 

58. 800 

4.713 

3. 959 

75 

63. 000 

3. 738 

3. 140 

80 

67. 200 

2. 759 

2.318 

85 

71. 400 

1.817 

1.526 

90 

75.600 

.982 

.825 

95 

79. 800 

.340 

.286 

100 

84. 000 

0 

0 

L. E. 

Radius: 1.704 Percent c= 1.431 inches 


Local unfairness of the model resulted in some minor varia- 
tions in the pressure distribution, as shown in figure 3, but 
these were not of sufficient magnitude to induce transition to 
turbulence. 

METHOD OF HEATING 

The heated portion of the wing was divided into four com- 
partments by spanwise bulkheads on which the heating ele- 
ments, nine rows of ordinary 120-volt incandescent lamp 
bulbs, were mounted (fig. 2). The inside surface of the 
aluminum skin was painted a dull black to increase the ab- 
sorption of radiant heat, but all other metal surfaces were 
bright. 

The heat input appropriate to each compartment to give 
a temperature difference across the boundary layer of 100° 
F. at a Reynolds number of 13,000,000 was calculated by the 
method of this report. The size, number, and location of 
the bulbs within each compartment were such as to give the 
most nearly uniform skin temperature possible within the 
practical limitations of the design. 

Power was supplied by a direct-current generator equipped 
with a remote voltage control which permitted a convenient 
means of adjusting the over-all applied voltage. 

Free-stream air temperatures were calculated from aver- 
age readings of three resistance-type thermometers located 
in the return passage a short distance ahead of the entrance 
cone. Adiabatic expansion through the entrance cone was 
assumed. 
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COMPUTATION OF HEAT TRANSFER 

METHOD 

The detailed analysis given in the appendix develops the 
following formulas, by which the local rate of heat transfer 
into both turbulent and laminar boundary layers may be 
computed. This method is applicable either to an airfoil 
or to a streamline body. The local rate for laminar flow is 

&r=0-700 ^ (t p — 1 0 ) (la) 

or the local heat-transfer coefficient is 

^=0.700 £ (lb) 

and for the turbulent flow, 


g*— 0.760 j- (t p —t 0 ) 
o T 


(2a) 


or the local heat -transfer coefficient is 

/>,= 0.760 A 

o r 


(2b) 


Heat transfer from an airfoil. — For the laminar boundary 
layer of an airfoil, 8 L is computed as in reference 2 from the 
pressure distribution as follows: 


3 / 17 p /c / u\ 8A7 , 

'c\uj Jo (v) d c (3) 


For the turbulent boundary layer of an airfoil, 8 T is com- 
puted as 



( 4 ) 


where f is the value of the turbulent boundary-layer para- 
meter as determined by a step-by-step solution of the rela- 
tionship of reference 3, given as 


di . 6.13 dU U 


dx' U dx 


=ir/(f) 


( 5 ) 



The value of / (f) is given in table II as taken from refer- 
ence 3, and may be plotted on semilogarithmic paper for 
ease in using. 


TABLE II 

NUMERICAL VALUES OF / (f) 

/(f) = 10.411 f-2 e-0.39Hf 


f 

10«/(f) 

f 

10«/ (f) 

14 

221.6 

24 

1.502 

15 

130. 39 

25 

.937 

16 

77. 47 

26 

.585 

17 

46. 11 

27 

.367 

18 

27. 97 

28 

.2308 

19 

16.96 

29 

. 1454 

20 

10.35 

30 

.0917 

21 

6.35 

31 

.0582 

22 

3. 906 

32 

. 0369 

23 

2.418 

33 

.0235 


7.S4.-147 


48 — 2 
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Heat transfer from a streamline body. — For the laminar 
boundary layer of a streamline body, b L is computed from 
reference 2 as 




( 6 ) 


The turbulent boundary-layer heat-transfer length may 
be computed as 

t'L 

(7) 


R, 


m 


where f is determined as for an airfoil from the equation 


dx 


6.13 dU 
U dx^ 


2.557 dr_U „ * 
r dx~ v m) 


(8) 


by the step-by-step process mentioned for calculating b T for 
an airfoil. 

Compressibility correction. — If the heat flow is to be ob- 
tained at free-stream Mach numbers such that the aero- 
dynamic temperature rise is an appreciable portion of the 
total temperature difference, a correction for aerodynamic 
heating should be made. The “heat-transfer temperature 
difference” to be used for a laminar boundary layer is 

(tp—to) = (t' — t 0 ) — 0.20M 0 2 T 0 [l-0.13 (jT)‘] (9) 

and for turbulent flow 

(t P - to) = (U -t 0 )-0 . 20 A/ 0 2 T 0 (10) 

where t e —t 0 is the desired temperature rise. 

The total rate of heat transfer from an airfoil or streamline 
body may be found as follows: 

1. Estimate the location of the transition point by the 
method of reference 4 for an increasing pressure gradient, 
or by reference 2 for a falling pressure gradient. 

2. Calculations for the laminar region ahead of the transi- 
tion point. 

(a) Compute the values of b L along the surface to the 
transition point by equation (3) for an airfoil, or by equation 
(6) for a streamline body. 

(b) With these values and the desired temperature distri- 
bution corrected for compressibility, compute the local rates 
of heat transfer along the surface by equation (la). 

3. Calculations for the turbulent region behind the transi- 
tion point. 

(a) From the value of b L at the transition point compute 
0, the momentum thickness as 


0 = 0.289 b L 

Using this value of 0, find the initial value of f at the 
transition point as 

f =2.557 log* 4.057(^) 

(b) With this initial value of f, calculate the values of 
f along the surface by equation (5) for an airfoil, or by 
equation (8) for a streamline body. With these values of 


f, compute b T along the surface by equation (4) for an 
airfoil, or by equation (7) for a streamline body. 

(c) Using these values of b T and the desired temperature 
difference across the boundary layer corrected for compres- 
sibility, compute the local rate of heat transfer along the 
surface by equation (2a). 

4. Integrate these local rates of heat transfer along the 
chord for both laminar and turbulent regions to obtain the 
total rate of heat transfer. 

HEAT-TRANSFER MEASUREMENTS 

As mentioned in the Introduction, a limited number of 
heat-transfer tests were made on the heated wing model to 
check the accuracy of the theoretical method. 

The experimental results are subject to several sources of 
error. All computed rates of heat transfer were based on the 
temperature distributions obtained at the center of the span, 
assuming that the spanwise variation was negligible. This 
was essentially true except for a small portion at each end of 
the wing. Precautions were taken to minimize the heat 
losses at the ends of the wing. These are not believed large 
since the design of the heating system allowed only slight 
transfer by convection, and the conduction of heat from the 
wing to its supports is negligible. Losses due to radiation 
from the wing have been computed as a maximum of 5 per- 
cent for the whole surface heated to 100° F. above the 
surroundings. This loss is not considered in the heat- 
transfer data. 

For the purpose of computing heat transfer, the chordwise 
temperature distribution was computed and plotted as 
heat-transfer temperature difference (t p —t 0 )', this is the ob- 
served temperature difference corrected for compressibility 
effect; that is, t v is the temperature measured by a thermo- 
couple in the skin minus the computed aerodynamic heating 
temperature rise. The value / 0 is the free-stream air temp- 
erature (distinguished from the local temperature just out- 
side the boundary layer of the wing, which will be higher or 
lower than t 0 due to adiabatic variations caused by the ve- 
locity field of the wing). 

RESULTS AND DISCUSSION 

Heat-transfer tests of the wing were made in two parts, the 
first concerned with comparing computed and measured 
values of the local rate of transfer from the airfoil surface to 
laminar and turbulent boundary layers, and the second 
with checking the total rate of heat transfer from the airfoil. 

The tests to measure the local rate of heat transfer were 
made at zero lift for two test Reynolds numbers, with free 
transition to obtain the heat-flow rate into a laminar bound- 
ary layer, and with transition fixed at 5-percent chord to 
determine the flow rate into a turbulent boundary layer. 
Figure 3 presents the pressure distribution over the wing at 
zero lift. The experimental procedure consisted in adjusting 
the heat input so that the skin temperatures were nearly 
constant along the chord. With this temperature distri- 
bution achieved, it was assumed that the second compart- 
ment of the wing, extending from 14.6-percent to 26.3- 
percent chord, was thermally isolated so that no flow of 
heat occurred in the skin or through the bulkheads. The 
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power input to this compartment was then measured by 
means of a voltmeter and an ammeter for comparison with 
the calculated rate of heat flow. 

Figure 4 shows that the desired constant chordwise temper- 
ature distribution was attained for the laminar boundary 
layer, but the data of figure 5 for turbulent flow indicate 
that while the distribution was nearly constant from 10- to 
30-percent chord, covering the region under consideration, 
the temperatures over the nose were excessively high. This 
came about through the heating difficulties resulting from the 
sudden change in heat-transfer coefficient at the point where 
transition was fixed. This type of distribution may have 
resulted in some change in the local values of the temperature 
gradient at the wing surface, though this effect should be 
small since no appreciable temperature gradient existed over 
the portion of the surface concerned, the boundries of which 
are indicated by the dotted lines in the figures. 

To obtain the computed values, the variation of the heat- 
transfer coefficient along the chord was calculated for each 
case as outlined under the section Method (results plotted in 
figs. 6 and 7). The heat input into the second compartment 



Figure 4 —Chordwise distribution of temperature difference, ( t p — to ), a , 0°; transition free. 



Figure 5— Chordwise distribution of temperature difference, ( t P — to ). a , 0°; transition 
fixed at 5 percent chord. 



Figure 6. — Chordwise distribution of heat transfer coefficient, h z . a , 0°; transition free. 



Figure 7.— Chordwise distribution of heat transfer coefficient, h x . a , 0°; transition fixed 

at 5 percent chord. 

was then computed from the heat-transfer coefficient and 
experimentally measured temperature difference (corrected 
for compressibility heating effects). Both the measured and 
the computed values of heat input into the second compart- 
ment are listed in the following table: 


(deg) 

RcX 10-« 

Boundary layer 

Calculated 
heat trans- 
fer (kw) 

Measured 
input (kw) 

Error 

(percent 

input) 

0 

6.71 

Laminar 

0.861 

0.875 

1.6 

0 

10. 72 

do 

1.048 

1.028 

1.9 

0 

6. 90 

Turbulent 

3.76 

3. 75 

.3 

0 

11.17 


4.01 

4.01 

.0 


The agreement is considered quite satisfactory, and is taken 
to indicate that the method for the computations of the heat 
transfer coefficient involves no serious errors despite the 
assumptions involved. 

Further tests for the purpose of establishing the validity 
of the method as regards the total rate of heat flow from a 
wing were made at a lift coefficient of 0.55 and 8.60 million 
Reynolds number. In addition to a test with free transition, 
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a second condition simulating; the formation of ice near the 
stagnation point by fixing transition at 5-percent chord on 
the lower surface was investigated. Heat input in both 
conditions was maintained at the maximum available from 
the apparatus, and no attempt was made to achieve a pre- 
determined temperature rise or cliordwise distribution. The 
chord wise temperature distribution obtained, corrected for 
the effects of compressibility heating, is plotted in figures 
8 and 9. 



Figure 8.— Chordwise distribution of temperature difference, (f p — <o). a, 0°; transition free. 



Figure 9.— Chordwise distribution of temperature difference, (t p —t o). a, 5°; transition 
fixed at 5 percent chord. 

The chordwise variation of heat-transfer coefficient was 
computed for both the transition-fixed and transition-free 
conditions, as outlined in Method. Tables III and IV 
present the computations for the transition-free condition 
and serve as an illustrative example. The pressure dis- 
tribution used for these calculations is given in figure 10. 
The value of heat input was then computed by use of the 
experimental temperature distributions (corrected for com- 


TABLE III 

COMPUTATION OF HEAT-TRANSFER COEFFICIENT 


Airfoil 65,2-016 c/=0.55 Rc=8.69X10« Mo-0.183 

Upper surface Transition at l^ c chord 0i% e = 154.5XlO- 6 fi% e = 18.25 


Position 

Ax 

(£) 

dU/V 

dx 

/(OX 106 

dj 

dx 

Af 

r 

10 6 5r 

h x 

i 

c 

(H) 

0.01 

0. 07 


1.61 





18 . 25 

166. 5 

0. 0191 



0. 07 


-0.667 

24. 6 

'si.'s 

3. 63 




.02 

. 14 


1. 57 





21.88 

245. 4 

’.'6136 



.07* 


-.572 

- 4 .V 7 

10. 36 

‘.‘73 




'.03 

". 2 i 


1.54 





22.61 

267.’ 1 

"imo 



. 07 ' 


’-.’443 

2. 92 

7.34 

’.’51' 




.04 

.28 


1.52 





23. 12 

283.2 

.0113 



'.07' 


-. 300 

’2.29' 

5. 53 

‘.’ 39 ’ 




.05 

".35 


1.50 





23.21 

296.9 

’.bios 



.14' 


-.270 


4.62 

. 65 




.07 

.49 


1.48 





24. 16 

31 7’ 4 

.bibb 



. 1 1 


-. 185 

"i.39~ 

3.32 

. If, 




.09 

.63 


1.44 





24.62 

335.3 

.’60952 



. 14 


*-.”l33 

T 12 ’ 

2. 57 

’ . 3(V 




. 11 

. 77 


1. 422 





24.98 

353. 4 

.‘66903 




.35 


-.069 

".’946 

" 1. 96 

.69 




. io 

i." 12 


1. 399 





25.’ 67 

379.5 

.‘66841 



. 35 


-.046 

".‘685 

1. 39 

‘.‘49’ 




.21 

1. 47 


1. 382 





26.16 

399.’ 6 

.00800 



.35 


-.0372 

".‘545 

i.’io 

.39' 




.26 

l . 82 


1. 386 





26.' 55 

4l6’6 

. 00778 



’.’35’ 


—.0286 

".’453 

’.’906 

’.’32' 




.31 

2.17 


1.356 






26. 87 

430.9 

.0074 i 



.35 


-.017i 

".'393 

’.‘739 

’.’26’ 




.36 

2.52 


1.350 





27. 13 

439. 8 

.00725 




. 35 


—.0143 

".’345 

.643 

’.’23’ 




.41 

2. 87 


1.345 





27.36 

448.’ 2 

.’66712 




.35 


-. 0343 

".'312 

‘.’677 

’.‘24’ 




.46 

'3.22 


1.330 




_ 


27.60 

461.0 

.’60682 




.35 


— . 09l5 

".’278 

.918 

‘.‘32' 




.51 

’3.57 

::::: 


:::::: 

:::::: 

::::: 

— 

27.92 

472.’ 1 

.’00676 


TABLE IV 

COMPUTATION OF HEAT-TRANSFER COEFFICIENT 


Airfoil 65,2-016 0=0.55 I? C =8.69X10« A/o=0.183 

Lower surface Transition at minimum pressure 


X 

c 

(*)* 

(?)“ 


(?)“ 

“*T 

10 4 $/, 

(ft) 

10* hx 

0.010 

0.0 







.020 

. 185 

0. 00100 

4X10— 6 

0. 00044 

0.745 

5. 22 

56.3 

.030 

. 290 

.0063 

32X10-6 

.0034 

.758 

5. 66 

52.0 

.040 

.370 

.0173 

144X10-6 

. 0105 

.862 

6. 04 

48.7 

.050 

.455 

. 0400 

400X10-6 

.0270 

.951 

6.66 

44.2 

.075 

.60 

. 124 

2460X10 -6 

.096 

1.225 

8. 75 

33.7 

. 100 

.71 

.247 

0. 0067 

.208 

1.401 

9. 80 

30.0 

. 150 

.85 

. 515 

.0254 

.475 

1.806 

12. 65 

23. 2 

.200 

.92 

.712 

. 0565 

.683 

2. 248 

15. 73 

18.7 

.250 

1.00 

1.00 

. 1006 

1.00 

2. 478 

17.35 

16.9 

.300 

1.07 

1.32 

. 1590 

1.36 

2. 670 

18. 70 

15.7 , 

.350 

1. 12 

1.59 

. 2310 

1.68 

2. 895 

20.27 

14.3 1 

. 400 

1. 16 

1.83 

. 3158 

1.97 

3. 13 

21. 90 

13.4 

.450 

1. 17 

1.90 

.4100 

2. 10 

3. 45 

24. 15 

12.2 

. 500 

1. 18 

1.97 

. .5080 

2. 14 

3.81 

26. 62 

11.0 


possibility) and the calculated heat-transfer coefficients of 
figures 11 and 12. The computed and measured heat input 
are compared in the table below. Results indicate satis- 
factory agreement. 


a (deg) 

PrX10-6 

Boundary layer 

Calculated 
heat trans- 
fer 

Measured 
input (kw) 

Error (per- 
cent input) 

5 

8. 69 

Upper surface 1 turbulent. 
Lower surface laminar. 

19.60 

23. 52 

16.7 

5 

8. 60 

Upper surface turbulent. 
Lower surface turbulent. 

22. 30 

23. 55 

5.3 


The experimental temperature distributions for the tests of 
the total rate of heat transfer show that the lieat-transfer 
temperature difference varies to a marked degree along the 
chord (figs. 8 and 9). This variation violates one of the 
assumptions underlying the development of the method; 
that is, that the temperature difference is constant along the 
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Figure 11 Chordwise distribution of heat transfer coefficient, h x , for upper surface, a, 
5°; transition free and transition fixed at 5 percent chord; R e , 8.69 and 8.G0X10 6 . 

chord, which must be true if the transfer of heat at all 
points along the surface is analogous to the transfer of mo- 
mentum. To what degree this assumption may be ignored 
has not been determined analytically, since the problem of 
considering the variation of temperature along the chord 

784347—48 



Figure 12. — Chordwise distribution of transfer coefficient, hx, for lower surface, a, 5°; 
transition free and transition fixed at 5 percent chord; Re, 8.69 and 8.60 X10 6 . 

presents difficulties which have so far prevented a solution. 
The experimental results, however, indicate that the ac- 
curacy with which the total rate of heat transfer can be 
computed is not greatly impaired by the temperature varia- 
tions experienced. Generalization of this result must await 
further experimental checks. 

The accuracy with which the local rate of heat transfer 
may be computed in a rising pressure gradient is dependent 
upon the accuracy with which the surface shear may be 
determined. Squire and Young’s method (reference 3) 
assumes that the turbulent boundary layer in a rising pres- 
sure gradient exhibits the same characteristics as the fully 
developed turbulent layer of a flat plate. The extent to 
which the relationship between the surface shear, the mo- 
mentum thickness, and the local velocity so derived remains 
valid is shown by the accuracy of the Squire and Young 
method in determining friction drag. It must be realized, 
however, that the method will fail if turbulent separation is 
imminent. 

The thickening of the turbulent boundary layer due to the 
rise in pressure acting on the displaced mass of fluid also is 
ignored by the assumption that the heat-transfer rate is 
proportional to the surface shear computed by Squire and 
Young’s method. Actually, the heat capacity of the bound- 
ary layer is increased by this thickening which tends to 
increase the rate of heat flow at the surface. This counter- 
acts the effect of the profile distortion, resulting from the 
same cause, which reduces the surface shear since it tends 
to cause separation. But this effect, too, is negligible for all 
cases where Squire and Young’s method may be applied. 

In concluding, it must be stated that while the method 
presented herein is subject to a number of broad assumptions 
in its development, the experimental evidence presented 
shows the total rate of heat flow may be calculated with 
reasonable accuracy. 
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CONCLUSIONS 

The accuracy of the method for determining the rate of 
heat transfer from an airfoil is shown to be good by the 
results of a limited experimental investigation. Since the 
correctness with which the heat transfer can be computed is 
dependent mainly on the accuracy with which the boundary- 
layer characteristics may be determined, it is expected that 
the method possesses the same accuracy for computing 
heat-transfer rates from a streamline body. 

Although the development of the heat-transfer formulas 
is based on the assumption that the skin temperature remains 
constant along the surface, the experimental results show 
that for moderate temperature variations the precision is 
still good. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 


APPENDIX 


DERIVATION OF THE METHOD 


1. Heat transfer into a laminar boundary layer. — -The 
theory of heat transfer into a laminar boundary layer was 
first investigated by E. Pohlhausen for the case of incompress- 
ible flow along a flat plate maintained at a constant temper- 
ature (reference 1). Pohlhausen’s solution is developed by 
solving the differential equation for the temperature bound- 
ary layer by using Blasius* solution for the velocity boundary 
layer. 

In order to arrive at a solution for an airfoil in an incom- 
pressible fluid, it is necessary 

(1) to assume that the temperature boundary-layer and 
the velocity boundary-layer profiles for the airfoil are related 
in the same manner as for Pohlhausen’s solution. (This is 
true if the temperature of the skin remains constant along 
the surface and if the thinning of the friction layer in a 
favorable pressure gradient due to the change in pressure 
acting on the displaced mass of fluid is negligible.) 


(2) to calculate the value of 
boundary layer and then determine 



for the velocity 
with the relation- 


ship resulting from (1). (The solution of the problem for 
the temperature of the skin varying along the chord has 
been prevented because the diffculties so far have been found 
insurmountable.) 

Pohlhausen’s expression for the temperature gradient in 
the boundary layer at the surface of the plate is given as 



The function a(a) is the first derivative of Pohlhausen’s 
function defining the temperature boundary layer which, for 


( 7 = 1 , is equivalent to the second derivative of Blasius’ func- 
tion for the velocity boundary layer. Pohlhausen found 
that a(a) is accurately given by the relationship 

a (a) = 0 . 664 ^ 

then 

( i),.,--°' 332 Wix 

Now, for the Blasius boundary-layer distribution 

(P) =0.332 J Up U 

\dy/v=o V 

so that 

(dt\ = _/^A y- (J P -t_o) 

\dy/y=o \ay/ y =o v U 

It is now necessary to determine for the airfoil at 

any chordwise position. For the Blasius profile 

/du\ 0.765 jj 
\dy)y=o~ h 

where 8 L is the thickness of the boundary layer where u= 
0.707 U. Substituting in the preceding equation 


(4£) =0.765 Vo- 

\dy/v - o h 

or, taking <r= 0.760 for air, the local rate of heat transfer is 


or the lieat-transfer coefficient is 

h x = 0.700 ~ 

The development in reference 5 of an expression for the 
heat-transfer rate based on Reynolds analogy gives results 
which are in complete agreement with the above if <7=1. 
However, the experimental results of reference 6 indicate that 
the expression \/0.760, that is, \/ a air , properly relates the 
velocity and temperature gradients in the laminar layer. The 
two methods give results within 10 percent of each other, 
which is sufficient for practical cases. 

The values of 8 L for laminar flow may be determined both 
for an airfoil and a streamline body by the method of reference 2. 

2. Turbulent boundary layer. The theory of heat transfer 
in eddying flow as given by Dryden (reference 7) requires 
the introduction of several new concepts. If the equations 
of motion for turbulent flow are written by placing 
u =u-\-u', v=v+v', w=w+w', where the bars indicate mean 
values and the primes indicate fluctuations, and these values 
substituted into the equations of motion for steady flow, 
similarity between equations so developed and the steady- 
flow equations can be shown by introducing a value of eddy 
viscosity e. Similarly, the concept _of eddy heat conduct- 
ivity, p, is introduced by placing t=t+t', in the equations 
of the temperature field. 
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These values of eddy viscosity, e, and eddy conductivity, 
P, however, do not have the same properties as p and k since 
they vary from point to point in the flow. Nevertheless, e 
and P can be shown to vary in the same manner from point 
to point in the fluid. This is done by introducing PrandtPs 
concept of a mixing length; that is, a length of path followed 
by a fluid particle before it becomes lost in the mass of 
eddying fluid. 

It is therefore shown if the shear, T ' = — pU 'y' } then 


or that 



da 

du 

dy 

dy 


e= pi 2 


du 

dy 


in which the mixing length l varies from point to point in 
the fluid. Now, if the eddy heat transfer is considered to 
be— c p v't', then the eddy heat conductivity is equal to 


Cppl 2 


du 

dy 


or equal to c p e, provided the mixing length for heat 


transfer is considered to be the same as the mixing length for 
the transfer of shearing stress. Dryden states that available 
experimental data show that the mixing lengths near a wall 
are closely equivalent for transfer of heat and momentum, but 
that the relationship falls down, for instance, in the wake of a 
heated body. Since the present case concerns heat transfer 
from a wall to eddying flow in a boundry layer, it is believed 
that this relationship is acceptable. 

For turbulent flow, it has been shown that the Prandtl 
number is equal to unity; that is, 


CpP- lurb ^ 

k turb 

If the Prandtl number is unity, then the thermal and dynamic 
boundary layers have the same profile (reference 6). If wc 
make the same assumptions as in step (1) for the laminar 
boundary layer, we may write 

(dt_\ = (ip — to) 

\dy) v= o \dy Jv-o U 


since the surface shear 


/du\ 
\dy) v-o 
(tp «o) 

— LpT 0 JJ 


This is the same formula developed by Reynolds for flow in 
pipes (reference 7). 

So it is seen that the problem of calculating the rate of 
heat transfer in turbulent flow is primarily a problem of 
calculating the surface shear along the airfoil. This may be 
done by the method of reference 3, in which Squire and 
Young write the relationship 


where 



f= 2.557 log e (4.075 (^) 


6 being the momentum thickness of the boundary layer. 
(This relationship is developed from von Iv&rman’s formula for 
the skin friction experience by a flat plate with a fully devel- 
oped turbulent boundary layer. This assumption becomes less 
and less true as the turbulent boundary-layer profile of an 
airfoil becomes changed in shape and approaches separation 
in a steep pressure recovery.) Substituting for r 0 


Qx — Cp 


P IP (t p- to) 

r 2 u 


or 


Qx Cp ^2 (tp ^ 0 ) 


kUpVc /4 

q x — c p p £ £-2 n y c (tp t' 0 ' 


Cj>jx k VcU 

k cr 2 v v 


(tp to) 


where p and k are values for uneddying flow. Substituting 
for C, j the value for air, 0.760, 


where U is the velocity outside the friction layer. (This 
relationship is dependent on the assumption that the tem- 
perature along the surface remains constant as for step (1) 
for the laminar boundary layer, and that the thickening of 
the boundary layer due to the increasing pressures acting 
on the displaced mass of fluid and the distortion of the profile 
thus resulting is negligible.) 


Now 



<Z 


r=0.760 7772 R c 


Ct 2 



^ 0 ) 


or considering 



a characteristic length for the 


turbulent boundary layer, 
then 

g*=0.760 (t v -t a ) 

0 j> 

and 

A r =0.760 ~ 

Of 

In calculating 8 T , f may be computed by the step-by-step 
solution of the equations of reference 3 for an airfoil or 
streamline body. 
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3. Compressibility effects on heat transfer. — Since the 
foregoing analysis has been made for incompressible flow, 
the effect of aerodynamic heating must be dealt with if the 
heat transfer is to be accurately obtained. The effect of 
compressibility may be considered, to a first approximation, 
simply as influencing the lieat-transfer temperature difference 
to be used in the above-developed equations; that is, a part 
of any desired increase in the skin temperature will result 
from aerodynamic heating, and this part of the temperature 
increase involves no expenditure of heat. 

The temperature field near the heated surface of an airfoil 
or body of revolution operating at high Mach numbers may 
be determined by superposing the heat-transfer temperature 
field on that due to the friction heating as in reference 8. 
Eckert (reference 9) has shown that the temperature field 
due to aerodynamic heating for cr= 1 may be expressed as 

<=fe+I=lW7i[l 

or 

MJT-^ M l *T l (JjJ 

Superposing this on the temperature field for heat trans- 
fer, which may be given as 

u 

t =z tp {tp to) ~JJ 

the combined temperature field is 

t=t p -(t p -t 0 ) (y)+ y T 2 1 MJT 0 M L >T L (£j 

It is evident that at y= 0 

(dt\ _ _ (l P _—t») / du\ 

\dy)y = o U \dyJv-o 

which indicates that the heat transfer corresponds to the 
heat-transfer temperature field, so that for compressible 
flow the only correction necessary is that of correcting the 
skin temperature for the rise in temperature due to aero- 
dynamic heating. 


Eckert has shown that the temperature of the surface due 
to compressibility effects t M , is 

V o' “ 2 M - T l 

tM=t 0 + T ^ Mo 2 r 0 -U^ m l *t l +j* 7 “ 1 m l *t l 

or 

t M =to+ y -^ m<?t 0 [l - (1 - v^) (^) 2 ] 

The heat-transfer temperature difference for the laminar 
region, <r= 0.760, is then 

(t p — t 0 ) = (t e — to) — 0.2 M 0 2 T 0 j^l 0.13 J 

and for the turbulent region, a= 1 

(tp-t 0 ) = (t-t 0 )-0.2M 0 2 T 0 
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Positive directions of axes and angles (forces and moments) are shown by arrows 


Axis 

Force 
(parallel 
to axis) 
symbol 

Moment about axis 

Angle 

Velocities 

Designation 

Sym- 

bol 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) 

Angular 

Longitudinal 

X 

X 

Rolling 

L 

Y >Z 

Roll.. 

<p 

u 

V 

Lateral 

Y 

Y 

Pitching 

M 

Z >X 

Pitch „ 

9 

v 

q 

Normal. ... 

z 

Z 

Yawing _ 

N 

X >Y 

Yaw 

* 

w 

r 











Absolute coefficients of moment 

p L p M 

Ll ~qbS ~%cS 

(rolling) (pitching) 


0 -2L 

° n ~qbS 
(yawing) 


Angle of set of control surface (relative to neutral 
position), 8. (Indicate surface by proper subscript.) 


4. PROPELLER SYMBOLS 


D 

Diameter 


V 

Geometric pitch 


p/D 

Pitch ratio 


V' 

Inflow velocity 


V, 

Slipstream velocity 


T 

Thrust, absolute coefficient 

n r 

CT -pn*D< 

Q 

Torque, absolute coefficient 

ri Q 

LQ ~ P n 2 D> 


P Power, absolute coefficient C P = 


c 6 

Speed-power coefficient 

<pV l 

Pn 2 

V 

Efficiency 


n 

Revolutions per second, rps 


4> 

Effective helix angle = tan -1 

Ga-m) 


1 hp =76.04 kg-m/s=550 ft-lb/sec 
1 metric horsepower =0.9863 hp 
1 mph = 0.4470 mps 
l mps=2.2369 mph 


5. NUMERICAL RELATIONS 

1 lb=0.4536 kg 
1 kg=2.2046 lb 
1 mi= 1,609.35 m=5,280 ft 
1 m= 3.2808 ft 



